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Abstract We have carried out CO J = 2 — 1 and CO J = 3 — 2 observations toward 
Tycho's supernova remnant (SNR) using the KOSMA 3m-Telescope. From these obser- 
vations we identified three molecular clouds (MCs) around the SNR. The small cloud in 
the southwest was discovered for the first time. In the north and east, two MCs (cloud A 
and cloud B) adjacent in space display a bow- shaped morphology, and have broad emis- 
sion lines, which provide some direct evidences of the SNR-MCs interaction. The MCs 
is revealed at -69^-59 km s~^, well coincident with the Tycho's SNR. The MCs asso- 
ciated with Tycho's SNR have a mass of ~ 2.13 x 10^ Mq. Position- velocity diagrams 
show the two clouds adjacent in velocity which means possible cloud-cloud collision 
in this region. The maximum value (0.66 ± 0.10) of integrated CO line intensity ratio 
(/co j=3-2/^co j=2-i) for the three MCs agrees well with the previous measurement of 
individual Galactic MCs, implying that the SNR shock just drove into the MCs. The two 
MCs have a line intensity ratio gradient. The distribution of the ratio appears to indicate 
that the shock propagates from the southwest to the northeast. 

Key words: ISM: individual objects (Tycho's Supernova Remnant (G120. 1+1 .4)) - ISM: 
molecules - supernova remnants 



1 INTRODUCTION 

When a supernova explodes near the molecular clouds (MCs), shock generated by supernova remnant 
(SNR) can accelerate, compress, heat or even fragment surrounding interstellar MCs. They also can 
enhance abundances with respect to quiescent cloud conditions of different molecular species. So the 
SNR-MCs interaction plays a very important role in the evolution of Interstellar Medium (Reynoso et al. 
2000). Moreover, the molecular lines observations of MCs adjacent to SNRs can shed light on the SNRs' 
dynamical evolution and physical properties. So far about half of the Galactic SNRs are expected to be 
in physical associated with MCs (Reynoso & Mangum 2001). The association is indicated by morpho- 
logical correspondence of molecular emission, presence of molecular line broadening within the extent 
of SNR, and presence of line emission with high high-to-low excitation line intensity ratio, etc (see 
Jiang et al. 2010). However, only few cases about the SNR-MCs interaction has been well confirmed. 
Even for some SNRs interacted with the surrounding MCs, the detailed distribution of environmental 
molecular gas is poorly known. 
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Tycho's SNR is known as a young type la SNR located in the Perseus arm (Albinson et al. 1986), 
with an age of about 400 years. The distance of Tycho's SNR is estimated to be 2.3 kpc (Kamper et al. 
1978; Albinson et al. 1986; Strom 1988). The remnant has been widely studied in various wavebands. 
From the neutral hydrogen HI observation, Schwarz et al. (1995) concluded that the distance should 
be 4.6=b0.5 kpc; Reynoso et al. (1999) have performed an HI absorption study toward the remnant, 
suggesting that Tycho's SNR is currently interacting with a dense HI concentration in the northeast 
(^L5i?=-51.5 km s~^), which locally slows down the expansion of the shock front. In radio continuum 
(Reynoso et al. 1997) and high-resolution X-ray wavebands (Hwang et al. 2002), they indicated that the 
radio continuum and X-ray emission along the northeastern edge of Tycho's SNR is strongest. Based 
on the VLA observation, Dickel et al. (1991) indicates that the northeastern edge of Tycho's SNR is 
expanding with a very small and decreasing velocity into the dense interstellar medium, and may interact 
with it. Lee & Koo (2004) made the CO J = 1 — line observation towards Tycho's SNR. They 
concluded that most of the CO J = 1 — emission around Tycho's SNR is between -67 and -60 km 
s~^, but the velocity component in intervals -63.5 ~ -61.5 km appears to be in contact with the 
Tycho's SNR along its the northeast boundary. Using CO J = 1 — line to study the environs of 
Tycho's SNR, Cai et al. (2009) found that the CO J = 1 — emission form a semi-closed molecular 
shell around the SNR, the emission in velocity range of -69 ~ -58 km is associated with the SNR. 
Because the calculated virial mass of clumps is greater than their gravitational mass, they suggested that 
CO molecular clumps are being violently disturbed by Tycho's SNR shock. 

In order to understand the evolution of SNR interacting with MCs and investigate the detailed 
distribution of the molecular gas around Tycho's SNR, we have performed CO J = 2 — 1 and CO 
J = 3 — 2 observation toward Tycho' SNR. The observations cover for the first time the whole area 
of Tycho's SNR in these frequencies. Due to the observed molecular lines at higher frequencies, we 
can attain higher angular resolution, which is critical to identify relatively compact core. Also, higher J 
transitions are relatively more sensitive to hot gas. Thus, we can detailedly understand the distribution 
of the molecular gas in the vicinity of Tycho's SNR. 

2 OBSERVATIONS 

The observations of Tycho' SNR were made in CO J = 2 — 1 and CO J = 3 — 2 lines using the 
3m KOSMA sub-millimeter telescope at Gornergrat, Switzerland in Jan, 2004. The half-power beam 
widths of the telescope at observing frequencies of 230.538GHz and 345.789GHz, are 130^' and 80'^ 
respectively. The pointing and tracking accuracy is better than 10". The system temperature were about 
120K. The medium and variable resolution acousto optical spectrometers have 1501 and 1601 channels, 
with total bandwidth of 248 MHz and 544 MHz, and equivalent velocity resolution of 0.21km 
and 0.29km s~^, respectively. The beam efficiency 5eff is 0.68 and 0.72 at 230 GHz and 345 GHz, 
respectively. The forward efficiency Fgff is 0.93. Line intensities were converted on the main beam 
scale, using the formula T^b = ^eff/^eff x ^a- 

Mapping observations are centered at RA(J2000)=00^25^22.7^ DEC(J2000)=64°08'50.44'' using 
the On-The-Fly mode, the total mapping area is 14.5' x 15' with a 0.5' x 0.5' grid. The observed 
data were reduced using the CLASS (Continuum and Line Analysis Single-Disk Software) and GREG 
(Grenoble Graphic) software. 

The 1.4GHz radio continuum emission data were obtained from the NRAO VLA Sky Survey 
(NVSS; Condon et al. 1998). 

3 RESULTS 

Fig.l (Left) shows CO J = 2 — 1 MCs distribution around Tycho's SNR. The SNR appear as complete 
shell on the whole and bright knots in southwest, which was detected by the 1 .4 GHz radio continuum 
emission. Three MCs were clearly identified around Tycho's SNR. Each has been designated alphabet- 
ically, cloud A, cloud B, and cloud C. In Fig.l (Left), cloud A-C are well associated with the 1.4 GHz 
radio continuum emission of Tycho's SNR. It appears that the cloud A and cloud B are spatially adja- 
cent. Here cloud C is for the first time identified around Tycho. Fig.l (Right) shows the CO J = 2 — 1 
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RA [J2000] Velocity (km s"^) 

Fig. 1 Left: CO J = 2 — 1 intensity map (gray scale) integrated from -69 to -54 km s~^, 
overlaid on the NVSS 1.4 GHz radio continuum emission contours (black contours). Gray 
scale levels are from 1.7 to 5.6 K km (the la value is 0.15 K km ). Two dashed lines 
mark the direction of the cuts in Fig. 4. Right: Spectra of molecular line CO J = 2 — 1 at the 
peak position of A-C cloud clumps. Letter A, B and C indicate the different MC clumps. The 
vertical lines in the spectra mark the peak velocity. 

spectra towards the peak position of Cloud A-C , respectively. From these spectra, we can see that there 
are several same velocity components in intervals -70 ~ -53 and -7 ~ 2 km s~^. The velocity component 
in interval -7 ~ 2 km should belong to the foreground source (Cai et al. 2009). A double-peaked 
feature is observed in intervals -70 ^ -53 for cloud A and cloud B, which can be fitted by Gaussian, re- 
spectively. Applying the rotation curve of Clemens (1985) together with Vq = 220 km s~^, where Vq is 
the circular rotation speed of the Galaxy, we obtained the distance of Cloud A-C. Tycho 's SNR exhibits 
a complete shell structure of about 8^ In view of the morphological correspondence of Tycho's SNR 
and the Cloud A-C, the tangent point (the center of the SNR) in the direction to the SNR is at ~ 4.9 kpc, 
which is rather close to the distance 4.6 ± 0.5 kpc obtained from the neutral HI absorption observation 
(Schwarz et al. 1995). The fitted and derived parameters for Cloud A-C are summarized in Table 1. In 
Table 1, the clouds have broad emission lines. For cloud A, the distance of component peaked at -56.8 
km is 4.4 kpc, which is smaller than that of Tycho's SNR (~ 4.9 kpc), so this velocity component 
can be related to foreground Galactic gas emission. Other components in intervals -70 ~ -53 will be 
analyzed below in relation with Tycho's SNR. 



Table 1 The spectral line parameters at the peak positions of clouds 



Name 


RA 


DEC 


Tmb (rms) 


FW (rms) 


Vlsr (rms) 


I (rms) 


Distance 




(h m s) 


(o / //) 


(K) 


(kms-^) 


(kms-^) 


(K kms"^) 


(kpc) 


Cloud A 


00 25 20.7 


64 13 20.44 


1.6(0.2) 


9.7(0.2) 


-63.7(0.1) 


6.7(0.2) 


4.9 




00 25 20.7 


64 13 20.44 


0.4(0.2) 


6.6(0.5) 


-56.8(0.2) 


1.0(0.2) 


4.4 


Cloud B 


00 25 38.7 


64 10 20.44 


1.1(0.1) 


9.7(0.3) 


-63.5(0.1) 


4.5(0.2) 


4.9 




00 25 38.7 


64 10 20.44 


1.2(0.1) 


6.3(0.1) 


-66.2(0.1) 


3.1(0.3) 


5.0 


Cloud C 


00 25 02.7 


64 06 20.44 


1.3(0.2) 


7.0(0.1) 


-63.6(0.1) 


3.0(0.1) 


4.9 



The present high- J transition observations enable us to map in detail the molecular gas toward 
Tycho's SNR. After a careful inspection of the CO component in intervals -70 ~ -53, we find that the 
only channel map can show the morphological and kinematical possible signatures of interaction be- 
tween Tycho's SNR and the surrounding MCs. Fig. 2 shows CO J = 2 — 1 intensity channel maps 
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Fig. 2 CO J = 2 — 1 channel maps each 1 km s~^, overlapping the NVSS 1.4 GHz radio 
continuum emission. The radio continuum emission is indicated as black contours. Central 
velocities are indicated in each image. 



over the velocity range of -70 ^ -53 km with interval 1 km s~^. The radio continuum emission is 
indicated as black contours. There is virtually no CO J = 2 — 1 emission at Vlsr > -54 km and 
^LSR < -69 km s~^, most of the emission in between -69 km and -59 km s~^. In addition, cloud A 
position has some faint emission in intervals -59 ^ -54 km s~^. The faint emission corresponds to a CO 
component peaked at -56.8 km s~^, as seen from the CO J = 2 — 1 spectra of cloud A in Fig. 1 (Right), 
which can be related to foreground Galactic gas emission from above analysis. Furthermore, cloud A 
position have some strong emission in intervals -67 ~ -59 km s~^, which can be related to a CO com- 
ponent peaked at -63.7 km s~^. This component showing an arc of molecular gas is striking coincident 
along the northeast boundary. In Fig. 2, there are two velocity components in the intervals -69 ^ -64 and 
-67 ~ -59 km for cloud B, corresponding to peak at -66.2 km and -63.5 km in spectra (see 
Fig.l), respectively. CO J = 2 — 1 component at -69 ~ -64 km appear to be in contact with the 
remnant along its eastern boundary. In addition, we can clearly see that the southwest CO J = 2 — 1 
molecular gas at -65 ~ -61 km corresponds to cloud C. 

The different transitions of CO trace different molecular environment. In order to obtain the inte- 
grated intensity ratio of CO J = 3 — 2 and J = 2 — 1 (/co j=3-2/^co j=2-i), we convolved the 8C of 
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RA [J2000] 

Fig. 3 CO J = 2 — 1 intensity map (dashed line) are superimposed on the line intensity ratio 
map (color scale), the line intensity ratios (/co 7=3-2/^00 j=2-i) range from 0.13 to 0.66 by 
0.06. The rms level is 0.10 (la). The wedge indicates the line intensity ratios scale. Black 
contours indicate the NVSS 1.4 GHz radio continuum emission. 



CO J = 3 - 2 data with an effective beam size of VlSO^ - 80^ = 102'' (Qin et al. 2008). The integrated 
intensities were calculated for CO J = 2 — 1 line in the same velocity range as for CO J = 3 — 2. The 
integrated range is from -69 to -54 km s~^. The rms noise level is 0.15 K km s~^(l<j). Fig. 3 shows the 
distribution of the ratio (color scale) overlaid with the distribution of the CO J = 2 — 1 line integrated 
intensity (black dashed contours); The radio continuum emission is indicated as black contours. The 
ratio values for the cloud A, cloud B, and cloud C are between 0.23 and 0.66. The rms level is 0.10 (la). 
The maximum value for the three MCs is 0.66, which is in agreement with the value measured in the 
Milky Way (0.55 ± 0.08; Sanders et al. 1993), in NGC 253 (0.5 ± 0.1 in the disk; Wall et al. 1991), and 
in M33 (0.69 ± 0.15; Wilson et al. 1997). It also appear near value measured in the starburst galaxies 
M82 (0.8 ± 0.2; Giisten et al. 1993). In Fig. 3, the high ratios (red color) are showed in the northeast part 
of cloud A and B. The ratio value in this region is greater that in the region where cloud overlaps with 
Tycho' SNR. The clouds have a ratio value gradient increasing from southwest to northeast, suggesting 
that SNR shock is expanding into cloud A and B, and start to compress the clouds. 

Assuming local thermodynamical equilibrium (LTE) for the gas and optically thick condition for 
CO J = 2 — 1 line, we use the relation N^^ ^ 10^A^co( Dickman 1978). The column density is 
estimated as (Garden et al. 1991) 

Nco = 6.9 X lO^^^^ii^:^ / T^^dv cm-2, (1) 



exp(-ll.l/Tex) 

The excitation temperature (Tex) of the CO J = 2 — 1 line is estimated following the equation Tex = 
ll.l/ln[l + l/(Tmb/ll-l + 0.02)], where Tmb is the corrected main beam temperature. Furthermore, 
their mass is given byM = jj^Nn^S/ (2.0 x 10^^), where the mean atomic weight of the gas /i is 1.36, 
S is the size of core region. The physical parameters of the core are summarized in Table 2. 

4 DISCUSSION 



Tycho 's SNR has complete radio shell on the whole and bright radio knots in southwest. Among the 
MCs surrounding Tycho 's SNR, we have seen that three clouds (cloud A, cloud B and cloud C) are 
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Table 2 The physical parameters of the core in LTE. 



Name 


R(pc) 


Tex(K) 


A^coCcm ^) 




M{Mq) 


Cloud A 


5.1 


5.7 


2.2x10^^ 


2.2x10^^ 


1151 




1.4 


3.8 


0.6x10^^ 


0.6x10^^ 


25 


Cloud B 


2.8 


5.0 


1.7x10^'^ 


1.7x10^^ 


275 




5.0 


5.1 


l.lxlO^'^ 


1.1x10^^ 


581 


Cloud C 


2.1 


5.3 


1.1x10^^ 


1.1x10^^ 


96 




-65 -60 
Vlsr (km s"^) 



-65 -60 
Vlsr (km s"') 



Fig. 4 P-V diagram constructed from the CO J = 2 — 1 transition for Cloud B. Left panel: Contour 
levels are from 0.3 K to 4.4 K by the 0.1 K, with a cut along the north-south direction. Right panel: 
Contour levels are the same as in left panel, with a cut along the east- west direction. 



spatially coincident with the SNR. Broad emission lines detected in cloud A and B and bow-shaped 
morphology suggest that the MCs may be interacting with the SNR. After a careful inspection of the 
CO component, we find that CO component of cloud A in intervals -67 ~ -59 km (peaked at -63.7 
km s~^) and CO component of cloud B in intervals -69 ~ -64 km (-66.2 km s~^) and -67 ~ -59 
km (-63.5 km s~^) are well associated with Tycho's SNR. The position- velocity diagram across the 
peak of cloud B along the north-south and the east-west direction, respectively, are constructed from 
the CO J = 2 — 1 line (see Fig.4). In Figure 4, obviously, the velocity components of cloud A and 
cloud B are adjacent. For cloud A, CO component in intervals -67 ^ -59 km may belong to cloud 
B, because cloud A and cloud B are adjacent in space and in velocity, and have the approximately same 
peak velocity in spectra. We concluded that cloud A and cloud B may not only simply be overlapping 
along the line of the sight, but also colliding each other. Furthermore, the new cloud C is also associated 
well with bright knots, the bright knots may be radio emissivity increase due to compression of the 
shocked material, suggesting that cloud C appears to be swept up by the SNR shock in the southwest 
area. 

The maximum value of integrated CO line intensity ratio for the three MCs is 0.66 ± 0.10. The 
comparison with previously published observations reveals that the /co j=3-2/^co j=2-i for the MCs 
associated with Tycho's SNR agrees well with the value measured in in the Milky Way (0.55 ± 0.08; 
Sanders et al. 1993), in NGC 253 (0.5 ± 0.1 in the disk; Wall et al. 1991), and in M33 (0.69 ± 0.15; 
Wilson et al. 1997). It also appear near value measured in the starburst galaxies M82 (0.8 ± 0.2; Giisten 
et al. 1993). The high /co j=3-2/^co j=2-i in starburst galaxies may be due to unusual conditions in 
these dense and hot regions (Aalto et al. 1997), while for normal MCs the most likely explanation is a 
significant contribution to the CO emission by low column density material (Wilson & Walker 1994). 
Moreover, the high /co j=3-2/^co j=2-i value (3.4) in the MCs interacting with SNR IC443 (Xu et 
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al. 2010) exceed previous measurement of individual Galactic MCs, implying that the SNR shock has 
driven into the MCs. For the MCs around Tycho's SNR, the /co j=3-2/^co j=2-i value (0.66 ± 0.10 
) may indicate that the Tycho's SNR shock just drove into the surrounding MCs. In addition, the MCs 
associated with HII regions have higher ratio of /co j=3-2/^co j=2-i and have higher gas temperature 
than those sources without HII regions, indicating that the high line ratio may be due to heating of the 
gas by the massive stars (Wilson et al. 1997). Different masers maybe occur in different astrophysical 
environment. The H2O masers are located near the MSX sources and within the maximum line intensity 
ratio /co j=3-2/^co j=2-i regions, suggesting that H2O masers occur in relatively warm environments 
(Qin et al. 2008). Thus, we concluded that the line intensity ratios based on the optically thick CO 
transitions indicate the temperature varies at different positions. The MCs associated with Tycho' SNR 
have a ratio value gradient increasing from southwest to northeast, may imply that a shock has driven 
into MCs and compressed MCs. Since the MCs are compressed, it lead to the temperature of molecular 
gas increases. Hence the line intensity ratio /co j=3-2/^co j=2-i will increase. We consider that high 
^co 7=3-2/^00 j=2-i is also identified as one good signature of SNR-MCs interaction system. From 
Table 1, the total mass of MCs is - 2.13 x 10^. 

5 SUMMARY 

We have presented large-area map of molecular clouds (MCs) in the vicinity of Tycho' SNR in CO 
J = 2 — 1 and CO J = 3 — 2 lines. The complete map suggests that MCs are mainly distributed along the 
northern and eastern boundary of the SNR. It also reveals a new cloud (cloud C) that appears to be swept 
by the SNR shock in the southwestern area. The bow-sharped morphology, the broad emission lines, as 
well as the integrated CO line intensity ratios (/co 7=3-2/^00 j=2-i) (0.66 ± 0.10) further suggest that 
these MCs are interacting with Tycho' SNR. The MCs is revealed at -69^-59 km s~^, well coincident 
with the Tycho's SNR. The MCs associated with Tycho's SNR have a mass of 2.13 x 10^ M©. In the 
northern and eastern MCs, we find a line intensity ratio gradient along the southwest-northeast direction, 
implying that Tycho' SNR shock has driven into the MCs. We concluded that the line intensity ratios 
based on the optically thick CO transitions indicate the temperature varies at different positions. Both the 
integrated intensity maps and position- velocity diagrams show the eastern and northern clouds adjacent 
in space and velocity. Together with the approximately same peak velocity in spectra, we found that the 
two MCs are interacting, indicating possible cloud-cloud collision in this region. 
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